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Abstract: We consider eigenvalue problems for general elliptic operators of ar- 
bitrary order subject to homogeneous boundary conditions on open subsets of 
the euclidean N-dimensional space. We prove stability results for the dependence 
of the eigenvalues upon variation of the mass density and we prove a maximum 
principle for extremum problems related to mass density perturbations which 
preserve the total mass. 
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£3 : 1 Introduction 



We consider a general class of elliptic partial differential operators 

0<H,|/3|<m 



^ ! subject to homogeneous boundary conditions on an open subset Q of R N with 



^ i mine measure, we assume max me coemcieins 



finite measure. We assume that the coefficients A n n are fixed bounded real- valued 



functions such that A a p = A$ a and such that Garding's inequality is satisfied. 
For such operators we consider the eigenvalue problem 

Cu = \pu , (1.1) 

where p is a positive function bounded away from zero and infinity. Problem 
(11. ip admits a divergent sequence of eigenvalues of finite multiplicity 

Ai[p] < •• < A n [p] < ... . 

In this paper we prove a few results concerning the dependence of A n [p] upon 
variation of p. 

Keeping in mind important problems involving harmonic and bi-harmonic op- 
erators in linear elasticity (see e.g., Courant and Hilbert |6J), we shall think of 
the weight p as the mass density of the body Q and we shall refer to the quantity 
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M = j n pdx as the total mass of Q. In the study of composite materials it is of 
interest to know whether it is possible to minimize or maximize the eigenvalues 
\ n [p] under the assumption that the total mass M is fixed (see e.g., Chanillo et 
al. [I], Cox and McLaughlin [7J [HI |9], Henrot [13]). In this paper we generalize 
the results proved in [TB] for the Dirichlet Laplacian. In particular, we prove the 
following maximum principle where we refer to non-zero eigenvalues: 

All simple eigenvalues and the symmetric functions of multiple eigenvalues of 
( HHP have no points of local maximum or minimum with respect to mass density 
perturbations preserving the total mass. 

See Theorem 14.31 for the precise statement. Moreover, we generalize a result 
of Cox and McLaughlin [8] and we prove that \ n [p] are weakly* continuous func- 
tions of p, see Theorem 13.11 This, combined with the above mentioned principle, 
implies that if C is a weakly* compact set of mass densities then for non-zero 
eigenvalues we have: 

All simple eigenvalues and the symmetric functions of multiple eigenvalues 
°f ( EUP admit points of maximum and minimum in C with mass constraint 
M = const and such points of maximum and minimum belong to dC. 

See Corollary 14.61 for the precise statement. The reason why we consider the 
symmetric functions of multiple eigenvalues and not the eigenvalues themselves is 
related to well-known bifurcation phenomena which prevent multiple eigenvalues 
from being differentiable functions of the parameters involved in the equation. 
Moreover, the symmetric functions of multiple eigenvalues appear to be natural 
objects in the study of extremum problems, see e.g., [T5J [TTJ HB]. In fact, in 
this paper we prove that all simple eigenvalues and the symmetric functions of 
multiple eigenvalues are real-analytic functions of p and we compute the appro- 
priate formulas for the Frechet differentials which we need for our argument, see 
Theorem 13.61 

Theorem 14. 31 and Corollary 14. 61 are proved for so-called intermediate boundary 
conditions in which case one of the boundary conditions is u = on dQ (see 
condition (14. ip and Example I2.12p . This includes the case of Dirichlet boundary 
conditions 

»=i = -=s^=°- on9sl < L2 > 

On the other hand, Theorems 13.11 and 13.61 are proved for a larger class of ho- 
mogeneous boundary conditions, including Neumann boundary conditions. See 
Remark 14.81 for a discussion concerning Neumann-type boundary conditions. 

Our work is inspired by known results by Krein [TJ], Friedland [TTJ, Cox and 
McLaughlin [H [8j [9] concerning the description of optimal mass densities sat- 
isfying the condition A < p < B in the case of the Dirichlet Laplacian. The 
expectation is that optimal mass densities are bang-bang solutions, i.e., minimiz- 
ers and maximizers satisfy the condition (p(x) — A)(p(x) — B) = on Q. Explicit 
solution to this problem was given by Krein [T3] for N — 1. Friedland [TTJ proved 
a general result for the minima of suitable functions of the eigenvalues on convex 
sets of mass densities. In Cox and McLaughlin [H [9] , among other results, it is 
proved that both the points of minimum and maximum of the first eigenvalue are 
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bang-bang. 

Our method is too general to give a precise description of the extrema. How- 
ever, our approach allows to state a maximum principle concerning all eigenvalues 
of a quite general class of elliptic operators which can be applied to arbitrary sets 
C of mass densities. 



2 Preliminaries and notation 

Let Q be an open set in and m G N. By W m ' 2 (Q) we denote the Sobolev space 
of functions in L 2 (Q) with weak derivatives up to order m in L 2 (Q), endowed with 
its standard norm defined by 

i 

||u||w».a(n) = (lM||2 (n) + 5^ ||-D a u|lz»(n)J , (2.1) 

a|=m 

for all u G W m ' 2 (Q). By W™' 2 {VL) we denote the closure in W m ' 2 (Q) of the space 
of C°°-functions with compact support in Q. Let V(Q) be a closed subspace of 
W m ' 2 (tt) such that the embedding V(Q) C L 2 {Vt) is compact. Let A a/3 G L°°(fi) 
be such that A a/3 = A^ a for all a, (3 G N$ with \a\, \(3\ < m. By 1Z we denote 
the subset of L°°(fl) of those functions p G L°°(Q) such that essinf^p > 0. Let 
p G 1Z be fixed. We consider the following eigenvalue problem 




in the unknowns u G V(Q) (the eigenfunction) and A G K (the eigenvalue). Note 
that problem ( 12. 2 p is the weak-formulation of problem ( II. ip subject to suitable 
homogeneous boundary conditions. The choice of the space V(O) is related to the 
boundary conditions in the classical formulation of the problem. For example, 
if V(fi) = W™ ,2 (Q) we obtain Dirichlet boundary conditions as in ( II. 2p . If 
V(Q) = W m ' 2 (tt) we obtain Neumann boundary conditions. UV(Q) = W m > 2 (tt)n 
W^ 2 {Q), for some k < m, we obtain intermediate boundary conditions. See 
Example 12. 121 below. See also Necas [T51 Chp.l]. 

It is convenient to the denote the left-hand side of equation (I2.2p by Q[u, <p). 
It is also convenient to denote by L 2 (Q) the space L 2 (Q) endowed with the scalar 
product defined by 

<Ui,u 2 > p = / UiU 2 pdx, V iti, u 2 G L 2 (VL). 
Jq 

Note that the corresponding norm ||tt|| z, 2 (0) is equivalent to the standard norm. 

We assume that the space V(fl) and the coefficients A a p are such that Gar- 
ding's inequality holds, i.e., we assume that there exist a, b > such that 

a \\ u \\w^Hti) < Q\ u i u \ + & lkllz,2 ( n) , (2-3) 

for all u G V(Q). Actually, in many cases it will be more convenient to normalize 
the constants a, b > in such a way that 

a \\u\\w^2 (n) < Q[u,u] + b\\u\\ 2 L 2 (n) , (2.4) 
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for all u G V(fl). For classical conditions on the coefficients A a p ensuring the 
validity of (12. 3p in the case of Dirichlet boundary conditions we refer to Agmon 
[U Thm. 7.6]. Moreover, we assume that there exists c > such that 

Q[u,u}<c\\u\\ 2 wm ,2 {n) , (2.5) 

for all u G V(fi). Note that since the coefficients A a p are bounded, inequality (12. 5 j) 
is always satisfied if is a bounded open set with Lipschitz boundary (actually, 
it is sufficient that fHs a bounded open set with a quasi-resolved boundary, see 
Burenkov [3J Thm. 6, p. 160]). 

Under assumptions (12. 4p . (12. 5p . it is easy to prove that problem (12. 2p has a 
divergent sequence of eigenvalues bounded below by —b. To do so, we consider 
the bounded linear operator L from V(Q) to its dual V(Q)' which takes any 
u G V(Q) to the functional L[u] defined by L[w][y>] = Q[u, <p\, for all ip G V(fl). 
Moreover, we consider the bounded linear operator I p from L 2 (Q) to V(fl)' which 
takes any u G L 2 (VL) to the functional 7 p [u] defined by =< u, ip > p , for all 

ip G V{VL). By inequalities (12. 4p . (12. 5 p and by the boundedness of the coefficients 
A a 0, it follows that the quadratic form defined by the right-hand side of (12. 4 j) 
induces in V(f2) a norm equivalent to the standard norm (12. ip . Hence by the 
Riesz Theorem, it follows that the operator L + bl p is a linear homeomorphism 
from V(fl) onto V(fl)'. Thus, equation (12. 2 j) is equivalent to the equation 

(L + 6J p )(- 1 )o/ p [u]=^ M (2.6) 

where 

/i=(A + 6)- 1 . (2.7) 
Thus, it is natural to consider the operator T p from L 2 (Q) to itself defined by 

T p :=io(L + bI p )^oI p , 

where i is the embedding of V(Q) into L 2 (Q). In the sequel, we shall omit i and 
we shall simply write T p = (L + bI P Y^ o I p . Note that 

< T pUl ,u 2 > p = I p [u 2 )[(L + bl p )^ o I p [ Ul \] 

= {L + bI p )[(L + bl p )^ o I p [u x ]][{L + fo/,)^ 1 ) o J p [ U2 ]], (2.8) 

for all u±, U2 G L 2 (Q). Thus, since the operator L+bI p is symmetric it follows that 
T p is a self-adjoint operator in L 2 (Q). Moreover, if the embedding V(Q) C L 2 (Q) 
is compact then the operator T p is compact. By inequality (12. 4p . T p is injective. 
It follows that the spectrum of T p is discrete and consists of a sequence of positive 
eigenvalues of finite multiplicity converging to zero. Then by (12. 7p and standard 
spectral theory, we easily deduce the validity of the following 

Lemma 2.9 Let Q be an open set in K Ar , m G N. Let A a p G L°°(Q) be such that 
A a p = Ap a for all a, (3 G Nq with \a\, \(3\ < m. Let V(fl) be a closed subspace 
of W m ' 2 (Q) such that the embedding V(Q) C L 2 (Q) is compact. Let p G 1Z. 
Assume that inequalities ( |£.^| ) and h2.5\) are satisfied for some a,b,c> 0. Then 
the eigenvalues of equation \2. Ip have finite multiplicity and can be represented 
by means of a divergent sequence X n [p\, n G N as follows 

fo Ei«i \B\<m A a aD a uDPudx 

X n [p}= min max Jn ^ a ^ m \ . (2.10) 

Ecv(n) u&e J n u 2 pdx 
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Each eigenvalue is repeated according to its multiplicity and 



K[p] >~b + 



a 



(2.11) 



for all n G N. Moreover, the sequence fi n [p] = (b + A n [p]) 1 , n G N, represents 
all eigenvalues of the compact self-adjoint operator T p . 

Example 2.12 We consider the case of poly-harmonic operators. Let m G N. 
Let A a/3 = S a 0m\/a\ for all a, (3 G with \a\ = \(3\ = m, where 5 a/ 3 = 1 if a = (3 
and 8 a p = otherwise. Let k G N , < k < m and V(tt) = W m > 2 (tt) n W fc,2 (^). 
A^oie i/iai ( [i^.^p and \2.5\) are satisfied for any b > where a, c > are suitable 
constants possibly depending on b. Moreover, if k = m and the open set Q has 
finite Lebesgue measure then the embedding V(Q) C L 2 (Q) is compact. If < 
k < m and the open set Q is bounded and has a Lipschitz continuous boundary 
then the embedding V(O) C L 2 (Q) is compact (actually it is enough to assume 
that Q is a bounded open set with a quasi- continuous boundary, see Burenkov J3, 
Thm. 8, p. 169]). Under these assumptions all corresponding eigenvalues X n [p] 
are well-defined and non-negative. 

Note that if k = m then V(Q) = W™ ,2 (VL) and by integrating by parts one can 
easily realize that the the bilinear form Q[u, tp] can be written in the more familiar 
form 



C = (— A) m subject to the Dirichlet boundary conditions hi. 2) . Recall that the 
Dirichlet problem arises in the study of vibrating strings for N = 1 and m = 1, 
membranes for N = 2 and m = 1, and clamped plates for N = 2 and m = 2. 
In the general case k < m, the classic formulation of the eigenvalue problem 

is 



where Bj are uniquely defined 'complementing' boundary operators. See Needs [Wf 
for details. For N > 2, m = 2 and k = 1 we obtain the problem 



which is related to the study of a simply supported plate. Here K is the mean 
curvature of the boundary ofQ. See Gazzola, Grunau and Sweers FTE/ for further 
details. 

Finally, we note that if m = 2 and k = then V(Q) = W 2 ' 2 (Q) and problem 
112. ty) is the weak formulation of the Neumann problem for the biharmonic operator 





= o, Vj = o,...,fc-i, 

BjU = 0, V j = 1, . . . , m — k, 



in f2, 
on <9f2 
on d£l 




on dQ, 
I s = 0, on dn, 

flu ' 




in f2, 
on <9f2 
on d£l 



(2.13) 
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which arises in the study of a vibrating free plate. Here div^n is the tangential 
divergence and Pqq the orthogonal projector onto the tangent hyperplane to dQ. 
See also Chasman J±$. 



3 Continuity and analyticity 

By the min-max principle (12.101) it follows that X n [p] is a locally Lipschitz con- 
tinuous functions of p G 1Z. In fact, one can easily prove that 

iw i \ r ii ^ min{A n [pi], A n [p 2 ]} + 26 

A„ [pi] - Xn[p2\\ < —r — t — f — r Pi - Pi x«(n) , 

mmjess mt pi, ess mt p2j 

for all pi,p2 G 1Z satisfying ||pi — p2\\L°°{n) < minjessinf pi,essinf p 2 }. In fact 
A n [p] depends with continuity on p not only with respect to the strong topology 
of L°°(Q) but also with respect to the weak* topology, which is clearly more 
relevant in optimization problems. The following theorem was proved by Cox 
and McLaughlin [8] in the case of the Dirichlet Laplacian and mass densities 
uniformly bounded away from zero and infinity. The proof can be easily adapted 
to the general case. Moreover, it is possible to replace the uniform lower bound 
for p by a weaker assumption. 

Theorem 3.1 Let £1 be an open set in R N , m G N. Let A a/3 G be such 

that A a p = Apa for all a, (3 £ Nq with \a\,\(3\ < m. LetV(Q) be a closed subspace 
of W m ' 2 (il) such that the embedding V(fi) C L 2 (f2) is compact. Let C C 1Z be a 
bounded set. Assume that there exist a, b, c > such that inequalities p^P and 
A2.5\) are satisfied for all p G C. Then the functions from C to R which take any 
p G C to A n [p] are weakly* continuous for all n G N. 

Proof. Since C is bounded in L°°(Q), it suffices to prove that given p G C 
and a sequence pj G C, j G N such that pj ^* p as j — > oo then A n [pj] — > A n [p]. 
To do so, we first provqj that for each n G N there exists L n > such that 
A n [pj] < L n for all j G N. Let n G N be fixed and u\, . . . ,u n G V(Q) be lin- 
early independent eigenfunctions associated with the eigenvalues Ai[p], . . . , A n [p], 
normalized by < u r , u s > p = S rs for all r, s — 1, . . . , n. Note that 



lim / u r u s pjdx = / u T u s pdo, 
Jn Jq 



for all r, s = 1, . . . , n. Thus 



n / 71 \2 p s 71 \2 

lim / I \ ^ r u r I pjdx — [ / lrU r I pdx, (3.2) 

uniformly with respect to 7 = (71, ... , 7„) G W 1 with |7| < 1. Let E be the linear 
space generated by u\, . . . , u n . By (13. 2 j) it follows that for any e > there exists 
j e G N such that 

In Ejam^rn K^ a uD^udx J n E|q|,|/3|< m A a pD a uD^udx 

L u2 Pi dx ~ L u2 P dx 

+e(A n [p] + 26) < X n [p] + e(A n [p] + 26) (3.3) 



1 This is clearly trivial if we assume that < a < p for all p £ C, in which case X n [p] < X n [a]. 
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for all u G E, j > j e . By combining (I2.10p and (13. 3p we deduce that X n [Pj] < 
\ n [p] + e(A n [p] + 26) for all j > j t , which implies the existence of a uniform 
bound L n as claimed above. The rest of the proof follows the lines of Cox [8]. 
Let u n [pj], n 6 N be a sequence of eigenf unctions associated with the eigenval- 
ues X n [Pj] normalized by < u n [pj],ui[pj] > p = 8 n i for all n, I G N. Note that 
Q[w n [pj], M n [pj]] = X n [Pj] for all j G N. By inequality (12. 4p . the sequence u n [pj], 
j G N is bounded in the space V(Q) equipped with the norm (12. ip . It follows that 
possibly passing to subsequences, there exists u n G V(VL) such that u n [pj] weakly 
converges to u n as j — > oo in V(Q), and there exists A n G K. such \ n [Pj] converges 
to A n as j — > oo. Moreover, since the embedding V(Q) C L 2 (Q) is compact we 
can directly assume that u n [pj] converges to u n strongly in L 2 (Q) as j —> oo. By 
passing to the limit in the weak equation 

Q[u n [pj], (f] = X n [pj) < u n [pj],(p > Pj , V (p G V(Sl) , 

it follows that A„ is an eigenvalue and of problem (12. 2 p and u n a corresponding 
eigenfunction. Note that < u n ,ui > p = S n i for all n, I G N, hence A n , n G N is a 
divergent sequence. It remains to prove that A n = X n [p\ for all n G N. To do so, 
assume by contradiction that there exists an eigenfunction u G V{Q) associated 
with an eigenvalue A of the weak problem (12. 2p such that < u, u n > p = for all 
n G N. Assume that u is normalized by ||m|| p = 1/(6 + A). By the Auchmuty 
principle [2J applied to the operator L + bl p , we have 

l Q[u,u] +b\\u\\ 2 L 2 (n) 

< 7, " \\ u ~ P u-i,pM\lI (O) , (3-4) 



2(6+ X n [ Pj \) 

for all u G V(Q) and n,j G N. Here P n ^ Pj u denotes the orthogonal projection 
in L 2 [p) of u onto the space generated by Wi[pj], • • • ,M„_i[pj] for all n > 2 and 
Po )Pj u = 0. By setting u = u and passing to the limit in (13 .4p as j — > oo, we 
obtain 



-276+X) " 2 — ~ m ^ ~ ~2{bTT) 

for all j G N, which contradicts the fact that A n — > oo as n — > oo. □ 

By classical results in perturbation theory, one can prove that X n {p\ depends 
real-analytically on p as long as p is such that X n [p] is a simple eigenvalue. This is 
no longer true if the multiplicity of X n [p] varies. In the case of multiple eigenvalues, 
analyticity can be proved for the symmetric functions of the eigenvalues. Namely, 
given a finite set of indexes F C N, we set 



Q[u,u} + b\\u\\ L2p{n) i 



K[F] = { P eTZ: \j\p\ ^ Xtlp], V j G F, I G N\F} 



and 



*fM= E X h[p\---K[pl h = l,...,\F\. (3.5) 

ji<—<3h 

Moreover, in order to compute formulas for the Frechet differentials, it is also 
convenient to set 

e[F] = {pe K[F] : X h [p] = X j2 [p], V jx,j 2 G F} . 
Then we have the following result 
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Theorem 3.6 Let £1 be an open set in M N , m G N. Let A a p G 6e suc/i 

t/icrf A a/ 3 = Ap a for all a, (3 <E Nq with \a\,\(3\ < m. Let V(ft) be a closed 
subspace ofW m ' 2 (fl) such that the embedding V(Cl) C L 2 (Q) is compact. Assume 
that there exist a, b, c > such that inequalities A2.3\) and \2.5\) are satisfied. 

Let F be a finite subset ofN. Then 1Z[F] is an open set in L°°(Q) and the 
functions A F ^ h are real- analytic in H[F]. Moreover, if F = U^ =1 F k and p G 
njJ =1 0[.Ffc] is such that for each k = 1, ...,n the eigenvalues Xj[p) assume the 
common value \F k [p) for all j G F^, then the differentials of the functions Ap,h at 
the point p are given by the formula 

n „ 

dk F)h [p][p\ = - J^c fe ^ / ufpdx, (3.7) 

k=l l£F k ^ Q 

for all p G L°°(Q), where 

0<hi<|Fi| V 7 j=l ^ 

0<h n <\F n \ 
h\-\ yh n =h 

and for each k = l,...,n, {ui}i eFk is an orthonormal basis in L 2 p (Q) of the 
eigenspace associated with A^Jp]. 

Proof. We set 

&F,h\p] = ( X n ip] + b)--- (A ih [p] + b) , 

ji,...,j h £F 
J Oh 

for all p G 7£[-F]. Note that by elementary combinatorics, we have 

LfM = E(-&)^('f L"/)^W , (3-8) 

where we have set A^o = A^o = 1. 

By adapting to the operator L + bl p the same argument used in [16] for the 
Dirichlet Laplacian, one can prove that 1Z[F] is an open set in L°°(0) and that 
Ai?/ l [p] depends real- analytically on p G Ti[F]. Thus, by ( 13. 8 p we deduce the 
real-analyticity of the functions Ap^. 

We now prove formula (13. 7p . First we assume that n — 1, hence F — F\ 
and p G Q[F]. For simplicity, we write Xf[p] rather than Ai?Jp]. The same 
computations used in [16] yields the following formula for the Frechet differential 
dA F>h [p] of A F)h at the point p G 

dA F , h [p][p] = -(X F [p] + b) h+1 ( 1 ^ ~ £ < rfT p [p]H, Mi > pj Vp G 

^ 7 leF 

(3.9) 
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By standard calculus and by recalling that T p U[ = (Xf[p) + b) 1 u i for all I G F, 
we have 



< dT p \p][ui],ui > p = -b<(L + hip) 1 dIp[p](L + bl p ) 1 I p u u ui > 



p 



+ < (L + blp) 1 dl p [p]u h ui > p = - A r F J P J - < (L + bl p ) l dl p [p\u u u l > f 

A F lp\ + 



(Mp}+b) 2 

hence 



ufpdx (3.10) 



n 



dk F]h [p][p\ = ~*f\p](Xf\p] +b) h - 1 ( lF h l _ f) £ / u iP dx ' ( 3 - U ) 



leF 

for all p G L°°(P). By (JUS) and (IXTTj) we get 

<*Aj?h[p][p] 

ft 



fc=l \ / \ / leF Jo. 

Xf[p] Ch-i) £ ( h k (Mp] + 6 )"(- & )^" fc E / 



which immediately implies f !3.Tj) for n = 1. We now consider the case n > 1. 
By means of a continuity argument, one can easily see that there exists an open 
neighborhood W of p in K[F] such that W C n£ =1 fc[F fe ]. Thus, 



0</ii<|Fi|,...,0</i„<|F n | fe=l 
/iiH \-h n =h 

on W. By differentiating equality (I3.12p at the point p and applying formula 
( 13. 7p for ra = 1 to each function A Fk ^ hk , we deduce the validity of formula ( 13 ,7p 
for arbitrary values of n G N. □ 

4 Maximum principle 

In this section we consider the case of general intermediate boundary conditions. 
This means that we assume that V(Q) is a closed subspace of W m,2 (Q) satisfying 
the inclusion 

V(Q) C W^\tt) . (4.1) 
Assume that Q has finite measure. For all M > we set 

L M = |p G L°°(tt) : j pdx = m| (4.2) 

The following theorem is a generalization of [THl Thm. 4.4] to the case of 
intermediate boundary conditions. 
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Theorem 4.3 Let all assumptions of Theorem 1 3. 61 hold. Assume in addition that 
Q has finite measure and inclusion Iji4-1\ ) holds. Then for all h = 1, . . . , \F\ the 
map Ap,h of1Z[F] R Lu to R which takes any p G H[F] R Lm to /ios no 

points of local maximum or minimum p such that Xj [p] that have the same sign 
and Xj[p] 7^ for all j G F. 

Proof. It is convenient to consider the real-valued function M defined on 
L°°(VL) by M[p] = J n pdx for all p G L°°(Q). Assume by contradiction the 
existence of p as in the statement. Then p is a critical point for the function 
Ap,h subject to the mass constraint M[p] = M. This implies the existence of a 
Lagrange multiplier which means that there exists c G R such that dAp^p] = 
cdM[p] (see e.g., Deimling [10, Thm. 26.1]). By formula (13. 7p . it follows that 



Jn \k=i ieF k J J a 




for all p G L°°(p,). Note that are non-zero real numbers of the same sign. Since 
p is arbitrary, it follows that 

= c, a.e. in Q. (4.4) 

Since ui G Wq' 2 (Q), then by a standard argument one can prove that the function 
(Sfc=i SieF fc ( A /l c fel' u ') 2 ) 1 ^ 2 belongs to the space W ' 2 (£l) and equals y/\c\ almost 
everywhere in f2. As is well-known the space W ' 2 (Q) does not contain constant 
functions apart from the function identically equal to zero. Thus c = and 
accordingly u\ = for all I G F, a contradiction. □ 

Remark 4.5 Theorem \4-3\ concerns mass densities p such that Xj[p] do not van- 
ish and have the same sign for all j G F . This assumption is clearly guaranteed for 
positively defined operators. Moreover, we note that the sign of the eigenvalues is 
preserved by small perturbations of p. Hence our assumption is not much restric- 
tive in the analysis of bifurcation phenomena associated with multiple eigenvalues 
different from zero. 

Finally, by Theorems 13.11 and 14.31 we deduce the following 

Corollary 4.6 Let all assumptions of Theorem \4-3\ hold. Let C C TV[F] be a 

weakly* compact set in L°°(Q). Assume that there exist a, b > such that in- 
equality {2.1$ is satisfied for all p G C. Let M > be such that C R Lm is not 
empty. Assume that the eigenvalues Xj[p] have the same sign and do not vanish 
for all j G F , p G C . Then for all h G {1, . . . , |F|} the map A F ^ h from C R L M to 
R which takes p G C R L M to A Fih [p] admits points of maximum and minimum 
and all such points belong to dC R L M . 

Proof. Recall that weakly* compact sets are bounded. Thus, by Theorem 13. II 
the functions Ap^ are weakly* continuous on C hence they admit both maximum 
and minimum on the weakly* compact subset C (ILm of C. By Corollary 14 . 3 1 the 
corresponding points of maximum and minimum cannot be interior points of C, 
hence they belong to dC R Lm- n 
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Example 4.7 Consider the poly-harmonic operators subject to Dirichlet or in- 
termediate boundary conditions as described in Example \2.12\ Let A,B& 
be functions satisfying the condition 

< ess inf A(x) < esssup£>(x) < oo. 

Let C = {p G L°°(Q) : A < p < B}. Clearly, C is a weakly* compact set. 
Moreover, since all mass densities p are uniformly bounded away from zero and 



infinity, inequality \2.1$ is satisfied for suitable constants a, b > not depending 
on p G C . Thus Corollary \4- 6] is applicable to all non-zero eigenvalues. It turns 
out that point of maximum and minimum p should coincide with A(x) or B(x) 
in a set of positive measure. 

Remark 4.8 Condition was used only to guarantee that V(Q) \ {0} does 

not contain constant functions. Thus, one may replace condition ^4- l\j by slightly 
more general conditions. For example one may assume that V{Q) C Wo'r(ft) 
where Wqp(Q) is the closure in W 1,2 (Q) of C°° -functions vanishing in an open 
neighborhood of a suitable subset of T of dVL. In this case, one would talk about 
mixed-intermediate boundary conditions. 

IfV(Q) is a closed subspace ofW m,2 (Q) containing constant functions differ- 
ent from zero, then we could argue as in the proof on Theorem \4-'J\ up to condition 
^4-4ty - Thus, in the general case one could simply characterize the critical mass 



densities of the functions Ap^ as those mass densities for which condition (j4-4\) 
is satisfied. Clearly, in the case of simple eigenvalues condition l[4-4\ ) reduces to 
u = const in Q which implies that A = 0. Thus, we conclude that the maximum 
principle stated in the introduction holds for all simple eigenvalues and all homo- 
geneous boundary conditions under consideration. As for multiple eigenvalues we 
note that the analysis of condition jjj4.4\ ) is not straightforward as it may appear 
at a first glance. Under suitable regularity assumptions on the eig en functions u\, 
u 2 associated with a double eigenvalue X of the Neumann Laplacian, one may 
prove that the condition u\ + u\ = const in Q implies that A = 0. However, we 
do not include such arguments here since we plan to perform a deeper analysis of 
Neumann and other boundary conditions in a forthcoming paper. 
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